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DETECTING AND MONITORING AGRICULTURAL VEGETATIVE WATER STRESS 


OVER LARGE AREAS USING LANDSAT DIGITAL DATA 
David R. Thompson and Oscar A. Wehmanen 


ABSTRACT 

A technique utilizing transformed Landsat digital data (for detection 
of agricultural vegetative water stress) was developed during the 1976 
South Dakota » U.S.A., drought. The procedure was expanded to the U.S. 
Great Plains during 1977 to evaluate the technii ue for detecting and 
monitoring vegetative water stress over large areas. The technique, Green 
Index Number (GIN), uses Landsat digital data from 3x6 nautical mile 
sampling frames (segments) to indicate when the vegetation within the 
segment is undergoing stress. At known growth stages for wheat, segments 
were classified as drought or non-drought areas. The remote sensing based 
information was compared to a weekly ground-based index (Crop Moisture 
Index) provided by the United States Department of Commerce. This 
cemparison demonstrated good agreement between the 18-day remote sensing 
technique and the ground-based weekly data. The moisture-stress detecting 
procedure developed over a small geographic area (South Dakota) was 
adapted to a larger geographic region (U.S. Great Plains). During the 
1977 crop year, Landsat data indicated that portions of the USSR and 
Australia wheat regions were undergoing drought. Drought isoline maps 
produced by GIN agreed with conventional sources. 
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INTRODUCTION 


Drought is a conmon occurence for many parts of the world and often 
adversely influences the affected nations* economy. Meteorological 
indexes %mich provide methods of computing the extent and severity of 
general drought have been developed. The Palmer drought index is the most 
satisfactory method for monitoring drought. During periods when a major 
drought is developing and spreading, the Palmer index combines 
precipitation and teuiperature data which provides a means for routinely 
assessing the areal distribution of the various degrees of drought 
severity (Palmer, 1968; Julian and Fritts, 1968). 

Palmer (1963) used a two-layer soil moisture model for estimating 
soil moisture deficiency on a month-by-month basis. These monthly 
variations are adequate for the study of drought; however, they do not 
reveal the shorter periods of soil moisture fluctuations which ar^ 
important during strategic phases of plant growth. The earlier drought 
study was modified by Palmer (1968) to calculate a weekly index of soil 
moisture requirements. This index uses the same two-layer mod^^ls to 
compute a weekly hydrological balance. It is published in the United 
States as the Crop Moisture Index (CMl) by USDA and NOAA. Although the 
CMI provides an indication of current plant-water relationships, it 
requires climatological analysis of a long record in order to derive the 
constants which define the moisture characteristics of the climate in the 
area of interest . 
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In many regions of Che irorld, Che necessary informaCion Co run Che 
CNI is noC available, noC reliable, or noC currenC. The use of remoce 
sensing for droughc moniCoring trould provide Che capability Co assess 
moisCure condicions anywhere in Che world on a regular basis in real cime. 
In Che Large Area Crop InvenCory ExperimenC (LACIE) (MacDonald, 

1976) , a Cechnique using LandsaC digical daCa Co deCecC rooisCure sCress 
was developed during Che 1976 SouCh DakoCa droughC (Thompson and Wehmanen, 

1977) . This Cechnique, Green Index Number (GIN), was evaluaCed in 1977 
over Che U.S. GreaC Plains and uCilized in real Cime in Che analysis of 
Che 1977 USSR and AusCralia wheaC crop. 

OBJECTIVES 

In 1976, using LandsaC digiCal daCa over a small geographic area, a 
procedure was developed for deCecCing and moniCoring agriculCural 
vegeCaCive waCer sCress. The objecCives of Chis sCudy were Co evaluaCe 
Che applicabil icy of Chis procedure over a large area and Co use Che 
Cechnique in operaCional moniCoring of Che 1977 USSR and AusCralian wheat 
crops. 
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APPROACH 


LACIE 5x6 nautical mile sampling frames (segments) were selected 
throughout the United States Great Plains (Figure 1). During the 1977 
crop year» these same segments were evaluated in real time and compared 
with identical segment data collected in 1976. The Green Index Numbers 
(GIN) were computed to indicate whether the vegetation was moisture 
stressed or normal. The 2-year remote sensing classification was compared 
to the CNI to evaluate the applicability of the technique over large 
areas. Also during the 1977 crop year, selected segments in USSR and 
Australia were monitored for stress. When stress wa3 indicated, all 
segmen«:8 in the region were processed and drought isoline maps were 
produced for these regions during the critical heading a^'d grain filling 
periods . 

Using ideas presented by Kauth and Thomas (1976), a screening number, 
the Green Index Number (GIN) was developed during 1976 drought in South 
Dakota (Thompson and Wehmanen, 1977). The Green Index Number is a value 
designated to sumnarize the condition of vegetation within a r.ampling 
frame. It is based on the ability to detect and measure the area of 
growing vegetation using all four Landsat bands and to measure its area 
within the region. In approximate terms, the GIN is the percentage of 
land in an area with a "healthy" cover of vegetation. The GIN is computed 
on a LACIE sample segment, an area (5x6 nautical miles) with 22,932 
Landsat image elements (pixels). The data is processed using an automated 
screening procedure which rejects pixels with values which are 
unreasonable for agricultural area. This procedure was established 
empirically after inspecting many LACIE segments. The procedure for 
computing GIN is defined by: ' 
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Landsat Channels 


*1 

*2 

For an observation where X ■ X 3 

X4 

Zl 

^2 

Z is computed, tdtere Z ■ Z 3 

Z4 

Z - RX 

where , 


.433 

.632 

.586 

.264 

-.290 

-.562 

.600 

.491 

-.824 

.533 

-.050 

.185 

.223 

.012 

-.543 

.809 


Each vector is inspected automatically, and any vector having 
values unreasonable for agricultural data is discarded using the following 
procedure: 

A pixel is accepted as good only if: 


Zj < 100 

-8 < Zj 

-19 < Zj 

-5 < Z, < 15 
4 
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Once the ecrcening has been performed, the histogram of the value 
truncated to integer, is accumulated for good pixels. Z^ is defined as 
the greenness channel in the Kauth-Thomas transformation and is a weighted 
difference of the infrared and visible channels. Bare soil has a low 
greenness level which changes with base level and sample segment location. 
The soil greenness, s^ for each segment is estimated to the greenness of 
that pixel greener than only 228 (approximately IZ) of the other good 
pixels. The greenness scale is normalised so that the estimated soil 
greenness is cero, thus, the green number, g, for a pixel is g ■ ~ s. 

The green number contains information about green vegetation. To compute 
GIN, the pixels with g ^ IS are counted, divided by 22,932 and multiplied 
by 100. The level 15 was observed empirically to represent healthy green 
agricultural vegetation. The GIN is an estimate of the percentage of 
pixels in a Landsat scene having green numbers high enough (>15) to 
indicate full cover of green vegetation. It is computed using only 
Landsat data. 
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It vaa detenained during the 1976 South Dakota drought that a plot of 
GIN veraua acquiaition tiae for a noraal, pradoaiinantly wheat aegaent 
ahould follow a predeterained curve (Figure 2). If an obaerved point for 
a aegaent fell into the ahaded region, the aegaent waa claaaified aa 
drought affected. The bounda for the ahaded region were defined 
eapirically with t defined aa the approxiaate apring eaergence in daya. 

For different areaa of yeara, the ahaded area ia aoved along the y axia 
from aide to aide to match the green-up curve. The bounda for the ahaded 
region are auch that a normal agricultural aegoient greena up at a rate 
greater than 1/2 percent per day for 40 daya and it exceeda a GIN value of 
20 for 30 daya. After thia time, the %iheat in a region haa completed 
grain filling and ia harveated and the greenneaa of the region may 
decreaae. Each acquiaition for a aegment waa plotted and the aegment waa 
claaaified aa moiature atreaaed or not moiature atreaaed. For compariaon, 
theae aegmenta were alao claaaified uaing the CMI for Crop Repotting 
Diatricta (CRD). A CRD waa claaaified aa drought affected if ita CHI fell 
below -0.5 for 2 conaecutive weeka. Both claaaif icationa were reatricted 
to similar time frames. It was poraible for a aegment to start normal, 
and then undergo moisture stress. If the GIN classification did not agree 
with the CMI at all times, the GIN was considered as not agreeing with the 
CMI (Table 1). 
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RESULTS 


U.S. Cre«t Plains Re«ulf 

The d«t« uted in Chi* iCudy consiaCed of 36 LACIE aegmenCa located 
throughout the U.S. Great Plain* wheat growing region (Figure 1). Theae 
aegiaenca contained at leaaC 5 per- <nt agricultural cropland. The final 
data aeC conaiated of 70 segment yeara^ for the 36 LACIE aegment* (Figure 
1, Table 1). The remote aenaing claaaif ication waa evaluated againat the 
Crop Moiature Index (CHI). Iltia index measure* the degree to which 
moiaCure requirement* uT growing crops were met during the previous week. 
The index is computed from average weekly values of temperature and 
precipitation. These values are used Co calculate, using previous soil 
moisture condition and current rainfall, Che actual moisture loss. If the 
potential moisture demand, or potential evapotranspiralion, exceeds 
available moisture supplies, actual evapotranspiration is reduced and Che 
CMI gives a negative value. However, if moisture meets or exceeds demand 
the index is positive. The CMI represents the average conditions over a 
several-county region (Crop Reporting District), so local moisture 
conditions may vary because of differences in rainfall distribution or 
soil types. The specific type of agriculture is not considered in the 
CHI, but it assumes a water-use curve typical of the leaf area index of 
the crops %fhich predominate in the region. 


A segment year is defined as an observation of one segment for a growing 
season. 
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The final data set of 70 segment years shows good agreement between 
the remote sensing classifications (GIN) and the CMI classification (Table 
1). The contingency table (Table 2), which compares the two 
classification methods, shows that the classification based on the CHI and 
GIN are the same 83X of the time. A chi-square test was employed to test 
the hypothesis that GIN and CHI are independent. The resulting test 
statistic, which under the above hypothesis foliows a chi-square 
distribution with one degree of freedom, was highly significant indicating 
that GIN and CHI are not independent. It was concluded that the GIN is 
detecting moisture through crop responses, and this procedure, developed 
over a small region, is extendable to larger areas. 

An inspection of the eleven disagreements on the classification 
results (Table 1) disclosr<* that the sell types at the segment location 
have a different water holding capacity than that used in the CHI model (5 
disagreements). Also, rainfall patterns at the segment location differ 
from the amount recorded at the weather stations used in computing the 
CHI. This is reflected in that the other 6 disagreements occutred in 
segments that were located on the edge of the CRD, and thus the CHI does 
not necessarily represent the conditions that exist at the segment 
1 ocation. 
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USSR AND AUSTRALIA RESULTS 


During the 1977 crop year, the Green Index Number (GIN) monitoring 
program was run selected I«ACIE 5x6 nautical mile sample segments 
located in the wheat regions of USSR and Australia. When moisture stress 
was indicated on these segments by the GIN programi all LACIE segments in 
these areas were monitored and drought isoline maps were produce* for 
Landsat passes during the critical grain filling stage. For the USSR| two 
passes during July represented this period, and for Australia, two passes 
during late September and early October. 

The GIN monitoring program showed that parts of the USSR spring wheat 
region was undergoing moisture stress during the Landsat passes of July 
2-18, 1977 (Figure 3). During this time, the GIN value should be 20 or 
greater for adequate moisture. The next Landsat overpass (Jul> 19* July 30, 
1977) showed that the area was still undergoing moisture stress (Figure 
4). A composite of the two GIN classification results show that much of 
the USSR spring wheat region experienced moisture stress during the most 
critical growth stage (Figure 5). While ground truth is not readily 
available for this area, the U.S. agricultural attache (personal 
communication) verified that this region had undergone moisture stress at 
that time. The official USSR small grains production estimates released 
were down considerably from normal (192 million metric tons vs. 225 
million metric tons) and thus also tend to verify that adverse conditions 
had occurred. 
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The critical growth stage for wheat in Australia ia during late 
September and October when the crop ia heading and the grain is beginning 
to fill. All lACIE segaients acquired for two Landsat overpasses during 
this tisM were run throu^ the CIN monitoring program. The Landsat pass 
for September 14-25 shows that smch of the area was undergoing mo{«i.jre 
stress (Figure 6). The next overpass (October 2-13) verified the area 
(Figure 7). A composite map shows that much of the wheat producing area 
of Australia was undergoing moisture stress (Figure 8). A drought review 
of Australia is produced by the Australian government each month. The 
October 1977 report shows the conditions that existed during the 
corresponding CIN interpretation of drought conditions (Bureau of Science, 
1977). These conditions are shown in Figure 9. The CIN classification 
represents only the area where LACIE segments are located. A comparison 
of Figures 8 and 9 shows general agreement. The Australia drought 
co'fivi.>ns are based on precipitation data and a compariaon with long-term 
averages. It does not take into account soil moisture available. An 
inspection of Landsat 100x100 nautical mile color IR images indicates that 
the vegetation is only marginally affected by the serious drought in the 
southeastern part of the country as reflected by the Australia drought 
report . 
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CONCLUSIONS 


The retui • of thle study demonstrate that a ''*hnique developed for 
sranitoring agricultural vegetative water stress over a small geographic 
region is expandable to large geographic areas and to different parts of 
the world. Thus, in areas of the world where ground truth is not 
available or reliable, it is possible to detect and determine the areal 
extent of stoisture stress using Landsat in an automatic sK>de. 
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TABLE 1. RESULTS OF GIN AND CMI CLASSIFICATIONS 


D ■ Moisture stressed conditions; W ■ Nonsal conditions; - ■ No data 


1977 1976 


Segment GIN 


1360 W 

1605 W 

1388 W 

1521 M 

1595 W 

1596 W 

1056 W 

1084 

1015 D 

1020 W 

1642 W 

1644 U 

1503 W 

1592 W 

1557 D 

1260 W 

1694 W 

1670 W 

1227 W 

1231 D 

1233 W 

1238 U 

1166 W 

1264 W 

1047 

1048 U 

1008 W 

1010 W 

1857 W 

1861 W 

1645 W 

1884 W 

1887 W 

1892 W 

1601 D 

1538 D 


CMI GIN CMI 

W D W 

WWW 
W - - 

W W D 

WWW 
WWW 

W D D 

W D D 

W D D 

W D D 

W W D 

W W D 

D D D 

WWW 
D D W 

W D D 

W D D 

W D D 

WWW 
W D W 

WWW 

WWW 
WWW 
W D D 

- D D 

W D D 

W D D 

W D D 

WWW 
W L W 

W W D 

WWW 
WWW 

WWW 
D W W 

D W W 
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TABLE 2. CONTIh *^.NCY TABLE OF GIN AND CMI CLASSIFICATION METHODS 


CMI 



Normal 

Stressed 

Total 

Normal 


5 

48 

Stressed 

6 

16 

22 

Total 

49 

21 

70 


X(,) - 27.89* 

*Sign if leant at the .001 level 
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8 1. Map of U.S. Great Plains showing locations of LACIE 5x6 

nautical mile sample segments. 
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Figure 3. Moisture conditions over USSR spring wheat from the Green 

Index Number (GIN) monitoring program (Landsat data acquired June 
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Index Nuaiber (GIN) monitoring program (Landsat data acquired July 




AUSTRALIA 





AUSTRALIA 



V 


u»dJ3 aq, Bcaj jpaq« ajipjjg,^ suoiaipuoa aanjsioH *8 a 



AUSTRALIA 





(JKJGINAi- ; 

Qfi EOGK QGAUa 




